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Abstract 

Motivated by recent experimental realizations of quintuple atomic layer films of Bi2Te3, the 
thermoelectric figure of merit, ZT, of the quintuple layer is calculated and found to increase by a 
factor of 10 {ZT = 7.2) compared to that of the bulk at room temperature. The large enhancement 
in ZT results from the change in the distribution of the valence band density of modes brought 
about by the quantum confinement in the thin film. The theoretical model uses ab initio electronic 
structure calculations (VASP) with full quantum-mechanical structure relaxation combined with a 
Landauer formalism for the linear-response transport coefficients. 
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In thermoelectric device applications a high value of the figure of merit (ZT) is desirable 

n ft 

for greater efficiency, [jj, |2| The thermoelectric figure of merit is defined as ZT = S 2 ctT/k, 
where T is the temperature, S is the Seebeck coefficient, a is the electrical conductivity, and 
k is the thermal conductivity which is the sum of the electronic (n e ) and the lattice (ni) 
contributions. For thermoelectric devices, a ZT value of above 3 or 4 at the room tempera- 
ture is required to be competitive with the conventional methods. {3] For the thermoelectric 
applications Bi2Te3 has become a material of particular interest, since it gives the highest 
ZT value of 0.68 at room temperature in its bulk form. |l| Recently, interest in this material 
las been further heightened by the fact that it is a three dimensional topological insulator 
4j, |5| , and it has been predicted that one- dimensional topologically protected modes in line 
dislocations can significantly increase ZT. 6] Many different approaches have been proposed 
and attempted to improve the thermoelectric performance of Bi2Te3, namely, by suppressing 



the lattice therma 



conductivity, [l-9] by tuning carrier concentration, by engineering the 



band structure, lOj and by reducing the device dimensionality. Ill, LL2[ The recent reports 
of the mechanical exfoliation and growth of quintuple layers of Bi2Te3 {13I \v\ and Bi2Se3 
15| motivate this study and underscore the importance of developing theoretical models to 
study the thermoelectric properties of few-atomic-layer thin films. 



The existing theoretical studies on the thermoelectric devices are most 



ther the continuum models 
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V based on ei- 



or the tight-binding methods. |19l . |20| Continuum 
models (effective mass approximation) though useful for their simplicity and computational 
efficiency, cannot capture the details of the atomistic effects that can be present in low- 
dimensional structures. On the other hand, tight-binding methods usually provide an ac- 
curate atomistic description of the bandstructures, however, these methods often employ a 
large number of empirical parameters obtained by fitting to experimental data or ab initio 
results that are available for the bulk materials. These parameters may not be transferable 
to low-dimensional structures where surface and interface issues can become important. In 
this study, we employ a state-of-the-art, density functional theory (DFT), ab initio method 
to calculate the electronic structures of bulk and two-dimensional (2D) Bi2Te3. The ther- 
moelectric parameters are then derived from the density of modes M(E) (i.e. a distribution 
of the conducting channels in energy) obtained from the electronic bandstructures. One 
advantage of the method is that it can be directly applied to low-dimensional structures. 
The calculations show that a room-temperature ZT value as high as 7.15 (around 10-fold 
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increase over the bulk value) can be achieved for atomically-thin 2D films of Bi2Te3. 

Fig. [1] shows the atomic structures of the two systems we study in this work: (a) a 
rhombohedral crystal structure of the bulk Bi 2 Te3 with the space group Dg d (R3m), it 
consists of five-atomic layers (Te-Bi-Te-Bi-Te) arranged along the ^-direction, known as 
quintuple layers (QL). The quintuple layers are separated from each other by weak van der 
Waals forces; (b) a free-standing film of Bi2Te3 with the thickness of 1 QL. In this work, all 
of the ab initio calculations (geometry optimizations and electronic structures) have been 
carried out using a projector augmented wave method within the framework of the Perdew- 
Burke-Ernzerhof (PBE)-type generalized gradient approximation of the density functional 
theory as implemented in the software package VASP. 2l| The Monkhorst-Pack scheme is 
adopted for the integration of Brillouin zone with a k mesh of 6 x 6 x 6 for the bulk and 
12 x 12 x 1 for the thin film. An energy cutoff of 400 eV is used in the plane wave basis. 
The optimized lattice parameters for the bulk are: a = 4.4842 A and c = 31.3767 A whereas 
the optimized thickness of the 1 QL thin film is d = 7.48 A. Spin-orbit coupling is included 
in the calculations of the electronic structure. 

Fig. [2] shows the ab initio electronic bandstructures of the bulk and QL film. The results 



23j | . The calculated 



are consistent with previous calculations for both bulk[22[ and thin films 
bandgap energy (E g ) of 90 meV for bulk Bi2Te3 is slightly smaller than the experimental 
bandgap of 150 meV. For the QL thin film, the bandgap increases to 190 meV due to quantum 
confinement effects. The underestimation of the bandgap energy in semiconductor materials 
with the DFT methods is a familiar issue. In order to achieve a better agreement with the 
experimental results we apply the so called 'scissors operator' and adjust the bandgap of the 
bulk Bi 2 Te 3 to match the experimental value by rigidly shifting the conduction bands up 
and the valence bands down by 30 meV around the midgap energy. The bandgap of the thin 
film is adjusted by the same amount. In this way, the calculated increase in the bandgap 
energy of the QL film compared to that of the bulk is left unchanged. After the rigid shift of 
the bands the bandgap energies become Eq = 150 meV for the bulk and Eg = 250 meV for 
the QL film. For the calculations of thermoelectric properties, we set the zero in the energy 
scale at the midgap energy for both the bulk and the QL film. 



The next step is to derive the thermoelectric (TE) parameters from the ab initio band- 
structures. The TE parameters are usually evaluated from the Boltzmann transport equa- 
tion (BTE). l] An alternative approach is the Landauer formalism 25, 26] which is more 
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convenient for performance comparisons between materials of different dimensions. 24J In 
this study, the objective is to compare the TE parameters between the 3D bulk and the 2D 
film, hence we have chosen the Landauer approach. Within the Landauer formalism in the 
linear response regime, the electronic conductivity (a), thermal conductivity for zero electric 
current (ft e ), and the Seebeck coefficient (S) are expressed as 



a = (2q 2 /h)I [1/fi-m], 
n e = {2Tk%/h){h - Il/h) [W/K - m] 



S 

with 



-(W?)£ [V/K], 



E — Ef 



T{E) 



dfo 



dE 



k B T J v ' \ dE 

where L is the device length and the transmission function T(E) = T(E)M(E) with M(E) 
as the density of modes (DOM). 26( In the diffusive limit, T(E) = X(E)/L with X(E) as 



the electron mean free path. 27| When phonon scattering is dominant, the mean free path 
can be written as X(E) = Xq, a constant. [20} The density of modes M(E) can be expressed 



as 
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27) 



M(E) = J2®(E-e(k ± )) (2) 

k± 

where G is the unit step function, and k± refers to all the k states in the first Brillouin zone 
perpendicular to the transport direction. Using Eq. (j2J), the DOM in any dimension can be 
numerically evaluated from a given e(k) simply by counting the bands that cross the energy 
of interest. Note that the expressions and the units indicated in Eq.fJTJ are same in all three 
dimensions provided the DOM is expressed in per unit area in each dimension. 

Fig. [3] shows the density of modes calculated from the ah initio bandstructures using 
Eq. (j2J). For the integration over the first Brillouin zone, k points are sampled on a uniform 
rectangular grid. The convergence of the final results has been ensured by using sufficient 
number of k points (101 x 51 x 10 for the bulk and 151 x 101 x 1 for the thin film). For 
the conduction band (E > 0), a slight increase in the value of the M(E) is observed for 
the QL film while the shape remains almost same as in the bulk. However, for the valence 
band (E < 0), the value of M(E) for the QL is much higher compared to that of the bulk. 
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Most importantly, at the valence band edge, the DOM for the QL turns on abruptly and 
has a peaked distribution. Previous work showed that a delta-shaped transport distribution 
function of the BTE (equivalent to M(E) in the Landauer formalism [^(J) maximized the 



ZT. 



16] 



Thermoelectric parameters for the bulk and the QL film are calculated at the room 
temperature from Eq.flT]) using the density of modes presented above. The results for the 
Seebeck coefficient (S) and the figure of merit (ZT) are shown in Fig. HI Values of the 
electron and hole mean free paths of Ao = 14 nm (for the conduction band) and Ao = 8 nm 
(for the valence band) give the best agreement with the bulk experimental data, and they are 
consistent with those used in previous Studi e,HQ FOTth e latti ee th e rmal conduct 
Ki we use the experimental value of 1.5 WK _1 m _1 . The results for the bulk Bi 2 Te3 show 
excellent agreement with the experimental data. For the thin film, there is a large increase 
in the thermoelectric parameters. The maximum in the Seebeck coefficient (S) increases by 
more than two times that of the bulk value, and the maximum ZT value of 7.15 obtained at 
Ep = —0.06 eV is around ten times higher than the bulk value of 0.68. The enhancements 
in the TE parameters are due to the improvements in the magnitude and the shape of the 
density of modes of the thin film brought about by the confinement in 2D. Note that for 
the thin film we use the same values for the parameters Ao and k\ as in the bulk. In reality, 
surface roughness, defects, and interface scattering present in 2D thin films can reduce those 
values. 

In conclusion, ab initio electronic structure calculations combined with a Landauer ap- 
proach for the linear-response transport coefficients show that the thermoelectric properties 
of the recently obtained Bi 2 Te3 quintuple layers can be significantly enhanced from those of 
the bulk. 
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Figure captions 

Fig. [TJ (color online) Atomic structures of Bi 2 Te 3 : (a) bulk unit cell; (b) free-standing slab 
with thickness of one quintuple layer (QL). Transport is assumed to be in the direction of 
the x axis (binary axis). 

Fig. [2j (color online) Ab initio bandstructures of Bi2Te3: (a) bulk; (b) thin film of 1 
QL thickness. The broken lines shows the position of the Fermi energy. For the bulk the 
symmetry points are given by: r(0,0,0), F(— |, 0, — |), L(0, |,0), and Z(|,|,|). For the 
thin film, the points are: T(0,0), A(~,|), and M(|,0). All the points are in the unit of 
{2ir/a) in the reciprocal space. 

Fig. EJ (color online) Density of modes M(E) per unit area for Bi 2 Te3: solid line for 
bulk and the broken line for 1 QL thick thin film. The midgap energy is set at E — 0. Note 
the sharp increase in M(E) for the thin film just below the midgap energy in the valence band. 

Fig. |H (color online) (a) The Seebeck coefficient (S) and (b) the figure of merit 
(ZT) at the room temperature (300K) as a function of the Fermi energy Ep for bulk and 
thin film of 1 QL thickness. The dotted black line shows the experimental data for bulk. 
The electron and hole mean free paths (Ao) are 14 nm and 8 nm, respectively, which are 
obtained by fitting our results to the bulk experimental data. The zero in the energy scale 
represents the midgap energy. 



8 




FIG. 1: 



9 




FIG. 2: 



10 




E(eV) 



FIG. 3: 



11 



(a) 
400 





-bulk 

-film (1 QL) ' 
-■ bulk Exp 


-0.2 -0.1 0.1 0.2 

E r (eV) 




-bulk 

—film (1 QL) , 
— bulk Exp 



-0.2 -0.1 0.1 0.2 

E r (eV) 



FIG. 4: 



12 



